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Recent results concerning the relation of topology and low-lying fermion modes are summarized. 



1. Introduction 

Low lying fermionic modes are believed to play 
an important role in QCD. This review covers 
some of the recent developments in the study of 
the relationship of low lying fermion modes with 
topology in QCD. The study of topology in gen- 
eral is not covered here. In particular, fermionic 
eigenmodes (including zero modes) are important 
in 3 areas discussed here: (1) low lying modes in 
particle spectrum, (2) low lying modes in global 
topology (e.g., chiral fermions), (3) low lying 
modes in the implementation of chiral fermions. 
The use of fermionic modes to probe for possible 
mechanism of chiral symmetry breaking in QCD 
are possible. The recent advent of the construc- 
tion of chiral fermions on the lattice are crucial 
in new studies. 



2. Spectrum 

2.1. Eigenmodes in Spectrum 

One expects low lying fermionic eigenmodes 
to be important in quark propagation because 
they should dominate the spectral sum for small 
enough quark masses. By computing some num- 
ber of the lowest lying eigenvectors of the Dirac 
operator one can expect that most of the physics 
of the quark propagator can be obtained. To 
make these statements more precise, an exam- 
ple is given of the computation of the rj' mass. 
Consider the spectral decomposition of a quark 
propagator - for simplicity use a hermitian form 
of the Dirac operator 



where D^^ is the standard Wilson-Dirac opera- 
tor, H is the hermitian Wilson-Dirac operator 
and the ^pi{x) are the eigenvectors with eigenvalue 
Xi. Then the spectral representation of the quark 
propagator is 



A,. 



(2) 



(1) 



The correlation function for the rj' involves a 
disconnected piece which is typically stochasti- 
cally estimated. A possible improvement is to 
truncate the spectral sum with the lowest few 
eigenvectors (which should give the largest con- 
tribution) and stochastically estimate the remain- 
der. The idea is to represent H — ^ ■ Hi + H± 
where H± are the remaining modes. Then in all 
relations where enters, find all the terms in- 
volving only Hi and the cross terms involving 
H± . Since the eigenvectors contain the full infor- 
mation (in the subspace) of the propagator (not 
just emanating from a single point), impose trans- 
lation symmetry by summing across the lattice 
for the lowest modes. This sum gives a volume 
times more statistics for the modes (hopefully) 
dominating the correlation function. For the re- 
maining terms involving H^^ use a stochastic es- 
timator. The linear system solutions involved in 
this latter step are accelerated since the condition 
number is lowered from projecting out the lowest 
eigenmodes. 

How many eigenmodes are needed? The scale 
should be set by the chiral condensate, hence the 
number should grow at least like the volume. 

Instead of computing the lowest eigenmodes, 
deflation can be used to compute the eigenvec- 
tors at the time of the linear system solution thus 
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Figure 1. Dependence of rj' effective mass on 
topological charge (Ref.||]). 
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Figure 2. One-parameter fit to a pure double- 
Goldstone pole form for the rj' hairpin correlator. 
The pion mass is fixed from the valence propaga- 
tor analysis (Ref.|jT^). 



accelerating the stochastic estimation 

One obvious question is whether to use 
for Wilson-like fermions the hermitian or non- 
hermitian operator. A comparison of the trun- 
cated spectral sum of the pion propagator shows 
poor convergence in the non-hermitian version 
when adding more eigenvectors The hermi- 
tian version shows slow but even convergence. For 
chiral fermions, the choice is irrelevant. 

The Wuppertal group has used the truncated 
spectral version of the hermitian quark propaga- 
tor for a calculation of the r]' correlation func- 
tion using Nf = 2 dynamical Wilson fermions |^] . 
A similar calculation was made for staggered 
fermions One goal is to determine the topo- 
logical charge dependence of the rj' mass. Us- 
ing a fermionic definition of the charge and bin- 
ning configurations according to Q, the effective 
masses show a Q dependence as seen in Fig. |l|. 
They find the truncated spectral sum method 
competitive with stochastic estimation for their 
large quark masses. The spectral sum method 
will not suffer dramatically going to smaller quark 
masses, however. The truncation of the spectral 
sum is evident in the not very cosh-like behav- 
ior of the r]' correlation function. The calcula- 
tion could be made exact by including the extra 
stochastic estimation. 



While not using a fermionic method, a calcu- 
lation by the UKQCD collaboration of the flavor 
singlet mesons with Nf — 2 dynamical 0{a) im- 
proved Wilson fermions has been made Q| . While 
statistically noisy, the OZI rule of suppressed sin- 
glet - non-singlet mass splittings is seen like in 
the quenched case 

3. Topology 

3.1. Topological Susceptibility 

The topological susceptibility is intimately tied 
with low lying fermion modes and is an important 
probe of the QCD vacuum relating to the Ua{^) 
problem. Pure gauge calculations of the topolog- 
ical susceptibility are consistent with the Witten- 
Veneziano prediction |6| . Another important test 
is to determine how the susceptibility depends on 
the quark mass. There are several recent large 
calculations of the topological susceptibility us- 
ing Nf = 2 dynamical fermions, from CPPACS 
(mean- field clover) Q , UKQCD (non-pcrtubative 
clover) g, SESAM/TxL (Wilson) §, and thin- 
linked staggered pO[ . By continuum arguments, 
the susceptibility is expected to behave like 

X(m) = ^(1 + 0{m)) = ^(1 + 0(m)). (3) 

Only the UKQCD calculation claims to un- 
ambiguously see the expected decrease of x("^) 
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Figure 3. Contributions to the ao propagator 
from an 77 — TT intermediate state. 



by lowering m. A comparison of the various re- 
sults at finite lattice spacing requires taking into 
account discretization errors p| . Naive linear 
extrapolations in (fixing F.^) give poor fits, 
suggesting that discretization effects are large. 
Besides large quark masses used there are other 
sources of systematic errors including the defini- 
tion of x(™) which involves contact terms and 
mixing with the unit operator For a chiral 
fermion action with a fcrmionic definition of 
these counter-terms are absent only in the Nf = 2 
chiral limit. They could be avoided by using a 
finite-volume scaling technique for small m 
applied to x(™)- 

3.2. Quenched Pathologies 

Another area low lying eigenmodes play an im- 
portant role are in quenched pathologies. How 
well is QCD described by an effective chiral the- 
ory of interacting particles (e.g., pions in chi- 
ral dynamics)? Suppressing the fermion determi- 
nant leads to well known pathologies as studied 
in chiral perturbation theory The patholo- 

gies are manifested in the rj' propagator missing 
vacuum contributions. A new dimensionful pa- 
rameter is induced in xPT. Power counting rules 
are changed leading to new chiral logs and pow- 
ers terms. The predicted pathologies for the pion 
mass were studied with Wilson fermions by CP- 
PACS . There are recent studies with Wilson 
fermions in the Modified Quenched Approxima- 
tion llqJTq] and using the overlap chiral fermion 
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Figure 4. Comparison of scalar ao propagator 
with the bubble sum formula fitted to the interval 
t=f-6 (Ref.@). 



operator 

The pion mass was one of the first observables 
studied for the effects of quenched chiral logs; 
however, it is not necessarily the easiest place to 
look. The additional quenched logs appear in the 
pion mass from the missing loop contributions in 
the r/' correlator. In ||l^, the unique piece of the 
77' correlator - the hairpin - was computed di- 
rectly and compared to quenched xPT. Fig. |^ 
shows the correlator along with the fit holding ttItt 
fixed to the simple mass insertion formula given 
by 

^ P g.P- (Tr75G(x, x)Tr75G(0, 0)) = 



(2^)4 
fp- 



-TO, 



2 



■ mt 



-fp 



(4) 



The correlator is well described by this form with 
the only free parameter the 77' mass insertion TOq. 
The pseudoscalar decay constant also shows the 
expected quenched artifact given by 



Jp 



quenched 



fp 



(5) 



while the axial vector decay constant does not 
show such a divergence. After considering vari- 
ous observables, the overall determination of (5 = 
0.065(f3) compared to S = 0.18 in Ref.||l|l. 

In Rcf. [|l^, the overlap-Dirac operator was 
computed on a 20^ lattice at a lattice spacing of 
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a — 0.13fm. Quenched chiral logs were observed 
in and fp giving 5 ~ 0.15 to 0.4. 

A dramatic example of quenching is exhibited 
in the non-singlet scalar particle For the ao, 
there is an ao ^ 77 — tt intermediate state as seen 
m Fig. I The vacuum contribution (a) is miss- 
ing in the quenched approximation along with the 
vacuum bubbles on the hairpin piece (c), but a 
remaining hairpin piece (b) is present. The re- 
maining hairpin piece has a quadratic divergence 
coming from the three pion propagators present. 
This quenched artifact causes a loss of positivity 
of the propagator as seen in Fig. |^. The fit is to 
the 1-loop bubble term as computed in quenched 
xPT and is in good agreement with the data jl6| . 
The difficulty in this Wilson fermion (MQA) ap- 
proach is the mixing of would-be zero modes and 
non-zero modes and in particular what are the 
contributions from each part. For example, is the 
lack of positivity a zero mode effect? It is not ex- 
pected to be, but more will be discussed on zero 
modes in the next section. 

While the ag quenched artifact can be viewed 
as the remnants of a decay, a true decay requires 
full QCD. The MILC collaboration has used an 
improved staggered fermion action in a Nj = 2+1 
calculation of the ao decay [|l^ . Evidence is found 
for a level splitting in the ao mass and an 77 -I- tt 
state. 

3.3. Condensate 

Several model calculations indicate the 
quenched chiral condensate diverges at T = 
(e.g., Ref. |19|). Damgaard |^^, shows via 
quenched xPT that the first finite volume correc- 
tion to the chiral condensate diverges logarithmi- 
cally in the 4-volume. Some simple relations for 
susceptibilities of pseudoscalar and scalar fields 
can be derived in the continuum and on the lat- 
tice with a chiral fermion action (see for example 
Rcf.[^). Define tt and ao fields by 

7r"(a;) = i'^{x)-i^T°-'il:{x),ao{x) = -i^{x)T''ip{x).{&) 

Then the tt and ao susceptibilities are given by 
1 



(^''(x)^"(0)) = -(VV') 

TO 

TO dm 



(7) 



Note in the quenched approximation the ag sus- 
ceptibility is given by the derivative of the chiral 
condensate. These relations are true including 
and excluding global topology terms. 

The Banks-Casher relation on a finite lattice is 



l^(7/;(x)^(x)) 



M 

mV 



(8) 



where / is some unspecified function of the non- 
zero eigenvalues. One sees the global topology 
term is irrelevant in the thermodynamic limit. 
Taking the appropriate limits the usual relation 
is recovered 



lim lim V(?/'(a;)V'(a;)) = 7rp(0+) 



(9) 



If the chiral condensate diverges without the 
topology contribution, the ao susceptibility must 
be negative and diverge. This divergence re- 
quires large enough physical volume to be ap- 
parent. Difficulties in the past include: (1) mix- 
ing of (would-be) zero and non-zero modes with 
staggered fermions, (2) mixing of topology and 
non-zero modes and contact terms with Wilson 
fermions, (3) until recently chiral fermion studies 
not on large enough lattices, e.g., random matrix 
model tests, spectrum tests, direct measurement 
tests. More direct tests of the divergence of the 
condensate will come from thermodynamics. 

3.4. Thermodynamics 

The deconfined phase of SU(2) quenched gauge 
theory, x 4, /3 = 2.4, above A'^t = 4 transition 
was studied using the overlap-Dirac operator Q . 
A study of the build-up of the density of eigenval- 
ues near zero, p{E) indicates the quenched chiral 
condensate (without topology) is diverging p^ . 
Define the density of zero eigenvalues p{E) from 
the derivative of the cumulative distribution 



N{E, V) = #(A > 0) where X<E 



p{E) 



dE V 



lim 



N{E,V) 



(10) 



As seen in Fig. |^, the cumulative distribution ap- 
pears to continually rise and track line on a log 
plot - hence the derivative (the condensate) di- 
verges with increasing lattice size. 
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Figures. Aplot of N{E,V)/{EV) from the small 
nonzero eigenvalues on four different lattices plot- 
ted on a log-log scale. The solid line is a least 
squares fit to the i = 18 data weighted by the 
statistical errors (Ref.[p2[). 



A diverging condensate indicates the spectral 
gap is closed. However, a decrease in the topolog- 
ical susceptibility is seen when crossing to T > 0. 
Has the quenched approximation distorted the 
vacuum so much that it has invalidated contin- 
uum arguments as to the nature of the T > 
phase transition or mechanisms of confinement? 
The claim is no - there is merely a large accu- 
mulation of near zero fermion modes. However, 
models predict changes in the vacuum structure 
crossing to the deconfined and (supposedly) chi- 
rally restored phase and low lying fermion modes 
can again be used to probe the vacuum. 

3.5. QCD Vacuum 

It is generally accepted that the QCD vacuum 
is characterized by strongly fluctuating gluon 
fields with clustered or lumpy distribution of 
topological charge and action density. Confine- 
ment mechanisms are typically ascribed to a dual- 
Meissner effect - condensation of singular gauge 
configurations such as monopoles or vortices (see 
review Ref. [^). Disagreements of various mod- 



els include: (1) instanton models provide chiral 
symmetry breaking, but not confinement, (2) cen- 
ter vortices provide confinement and chiral sym- 
metry breaking, (3) composite nature of instan- 
tons (linked by monopoles - calorons) at T > 0. A 
different approach relying on properties of gauge 
field correlators and not specific gauge field types 
describes a large class of non-perturbative phe- 
nomena |2^. 

Singular gauge fields are probably intrinsic to 
SU(3) (e.g., not just Gribov copies associated 
with gauge fixing). These singularities impose 
boundary conditions on quark and gluon fluc- 
tuations which moderates the QCD action. For 
example, instantons have locked chromo-electric 
and magnetic fields i?" = ±_B° that decrease in 
strength in a certain way. 

In a hot configuration one expects huge contri- 
butions to the action beyond such special type of 
field configurations. However, there could possi- 
bly be domains of (near) field locking. In recent 
calculations j25| , these domains have been shown 
to be sufficient to produce chiral symmetry break- 
ing, and confinement (area law). 

3.6. Instanton Dominance in QCD? 

Wittcn has argued |2^] that topological charge 
fluctuations are clearly involved in solving the 
[/yi(l) problem. However, the dynamics of the 
rj' mass generation need not be associated with 
semiclassical tunneling events. In particular, 
large vacuum fluctuations from confinement also 
produce topological fluctuations. Large Nc ar- 
guments are incompatible with instanton based 
phenomology - namely instantons produce and rj' 
mass that vanishes exponentially, but from large 
Nc chiral dynamics suggest that m^y ^ l/Nc- 
Witten speculated that the ij' mass comes from 
the coupHng of Ua{1) anomaly to topological 
charge fluctuations and not instantons. 

3.7. Local Chirality 

A test of Witten's conjecture was put for- 
ward with a local measure of chirality of non-zero 
modes ||2^. The relative orientation of the left 
and right handed components of low lying eigen- 
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Figure 6. Chirality histograms for the fowest two 
non-zero modes of the overlap-Dirac operator at 
the 2.5% sites with the largest tlj''ip{x) on the six 
ensembles with Iwasaki gauge action. The sys- 
tems in the left column all have approximately 
the same volume in physical units. The systems 
in the top two panels in the right column also 
have the same, about a factor of 3 larger, volume, 
while the system in the lower right hand corner 
panel has the largest volume (Ref. [p9[). 



modes is used 



tan(^-(l+X(x))) - 



(11) 



Here X{x) is a site dependent measure of the lo- 
cal chirality. A completely chiral state (an exact 
zero mode) have for some sites X{x) — ±1. These 
regions of local chirality come from (near) lock- 
ing of the chromo-electic and magnetic {E and B) 
fields. Wilson fermions were used as a test p7| , 
and since the chirality appeared random the claim 
was there is no instanton dominance in the QCD 
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Figure 7. The continuum topological susceptibil- 
ity in units of the string tension plotted against 
1 /N"^ . A linear extrapolation to iV = oo is shown 
(Ref. 0). 



vacuum. A flurry of papers using improved im- 
proved actions followed [ p8| , p9[ . 

The chirality histogram for the overlap-Dirac 
operator is shown in Fig. ^. Shown is the his- 
togram of X{x) for 2.5% of the sites with the 
largest ipip at three physical volumes keeping only 
the lowest two non-zero modes. Keeping more 
small non-zero modes indicates a finite density 
of such chiral peaked modes which may survive 
the continuum limit |^9|] . The mixing (trough) 
in the histograms are not related to disloca- 
tions. Namely, comparing histograms from dif- 
ferent gauge actions and vastly different gauge 
field dislocation content appear identical indicat- 
ing mixing comes from large sized eigenmodes. 
In addition, no significant peaking was found in 
a study of C/(l) in the confined phase where in- 
stantons are not present 

These results are all consistent with instanton 
phenomology; however, more generally they are 
consistent with suitable regions of (nearly) locked 
E & B fields. 

3.8. Large Nc 

Why is there a concern of reconciling the large 
Nc limit of QCD with studies of the QCD vac- 
uum? The large Nc description of QCD is quite 
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successful phenomenologically. It forms a ba- 
sis for the valence quark model and OZI rule, 
and well describes hadron spectra and matrix el- 
ements. The Witten-Veneziano prediction for re- 
lating the topological susceptibility and meson 
masses may well be reached at Nc = 3. 

How do gauge theories approach the limit? The 
prediction is that for a smooth limit, one should 
keep a constant t'Hooft coupling, A = g'^N as 
Nc —>■ oo. Is the limit realized quickly? Fig. |^ 
shows a test for the quenched gluonic topological 
susceptibility The large N^. limit is appar- 

ently realized quickly (seen more definitely in a 
2-1-1 study and is consistent with the pre- 

dicted scaling Q. 

Revisiting the issue of local chirality, it was 
found that the amount of peaking observed for 
example in Fig. ^ decreases at a fixed scale set 
by the string-tension as Nc increases from 3 to 
4 ||3^. This disappearance of local chirality is 
consistent with large (instanton-like) modes dis- 
appearing. In fact, Witten (26| predicts a strong 
exponential suppression of instanton number den- 
sity. Teper argued that including combina- 
torial factors for the measure could compensate 
for the exponential suppression leading to a sharp 
cutoff of the size distribution for small modes, but 
non-zero above the cutoff. However, given this 
limited step in increasing Nc the decrease in local 
chirality is consistent with the large Nc predic- 
tions. More work is needed to sort out contribu- 
tions to the size distribution from dislocations. In 
any case, the large Nc limit is an interesting place 
to gain additional understanding of the QCD vac- 
uum. Chiral fermions have been essential in these 
studies. 

3.9. Eigenmode Dominance 

How much are hadron correlators dominated by 
low eigenmodes? In recent work [|5[|6| compar- 
isons were made of the truncated and full spec- 
tral decomposition of various correlation func- 
tions using the (APE smeared) overlap-Dirac op- 
erator. The lowest 20 modes (including zero 
modes) were computed. These new results with a 
chiral fermion action draw some of the same con- 
clusions as older works HJs^. Namely, for small 
quark masses the non-singlet pseudoscalar corre- 



lator (pion) is well approximated by the truncated 
spectral sum. The non-singlet vector correlator 
only saturates at long time distances. These re- 
sults are consistent with instanton phenomology. 

The QCD sum rule approach parameterizes 
short distance correlators via the operator prod- 
uct expansion and long distance by condensates. 
There are large non-pertubative physics in non- 
singlet pseudoscalar and scalar channels. With 
the overlap-Dirac operator above |36), the trun- 
cated spectral sum for the point-point propaga- 
tors shows the appropriate attractive channel for 
the pseudoscalar and repulsive channel for the 
vector predicted by instanton phenomology and 
consistent with the Wilson case p7| . 

R,{x) = n,{x)/U'^{x), n,(x) - Tr(Jf (x)Jf (0)), 
Jf(a;) =i^(x)T"r(z)i^(a;) (12) 

Fig. |8| shows the ratio Ri{x) for the full and trun- 
cated pseudoscalar case with the free propagator 
for different quark masses. Saturation requires 
few modes for the lightest masses. One caveat 
is the use of APE smearing for the gauge links 
which could adversely affect the short distance 
part especially for heavier masses. 

3.10. Screening Correlators 

In the high temperature phase of QCD, one 
expects restoration of chiral symmetry exhibited 
in the limit ruqa — > by the equivalence of (non- 
singlet) screening correlators 

Cs{z)^~Cps{z), Cv{z)^Cav{z) (13) 

summed over the dimensions orthogonal to the 
z direction. A quenched calculation |Q with 
the overlap-Dirac operator at the A^t = 4 tran- 
sition shows good agreement in the vector (V) 
and axial-vector (AV) channels whereas previous 
Nf = & 2 calculations did not show agreement 
in the scalar (S) and pseudoscalar (PS) channels. 
The zero mode contributions tend to spoil Eq. |l^ 
but are a finite volume artifact. One can sub- 
tract them out directly or cancel them by com- 
paring {Cs{Z) — Cps[z))l2 to Cps{z) without ze- 
romodes. The agreement without the zero modes 
tends to support the conclusion that chiral sym- 
metry is restored and parity doubling is seen. 
However, the lack of a diverging chiral condensate 
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lap and domain wall (DWF) constructions see 
Refs. [41 S]. In particular, the five-dimensional 





Figure 8. Saturation of the point-to-point pseu- 
doscalar correlator by low-lying eigenmodes of the 
hermitian overlap operator. Results for 4 quark 
masses are shown with mps/mv ranging from 
0.34 in (a) to 0.64 in (d). Octagons show the full 
hadron correlator. Squares show the contribu- 
tion from the lowest 10 modes. Diamonds show 
the contribution from the zero modes, which just 
scales as l/{amqY (Rcf. pGf). 

seen in Rcf. probably from a small physical 
spatial volume indicates that increasing the size 
could break the parity doubling once again. More 
work to resolve this issue is needed. 

3.11. Localization 

With chiral fermions, one can test the issue of 
localization of fermionic eigenmodes and can pos- 
sibly compare to localization of gluonic quanti- 
ties. The latter is however strongly affected by 
dislocations in the gauge fields forcing the need 
for cooling or smearing. This unfortunately can 
affect part of what one was hoping to test. How- 
ever, fermionic studies p5| , p6|j3j - 40 provide evi- 
dence for localization. 

4. Chiral Fermions 

While low eigenmodes of chiral fermions have 
been discussed, the implementation of the chi- 
ral fermion operators is strongly affected by low 
lying modes. For recent reviews of the over- 



domain wall operator is equivalent to the four- 
dimensional overlap after taking the extent of the 
fifth dimensional to infinity. The overlap operator 
has the form £>o«eriap(0) = (1 + j5s{H{-M)))/2 
where H{—M) could be the super-critical hermi- 
tian Wilson-Dirac operator with mass M . The 
numerical implementation of e(i?(— M)) whether 
by an extra dimension via DWF or some ra- 
tional approximation is adversely affected by 
small eigenvalues of H{—M). The deviation 
of e{H{—AI)) from ±1 induces chiral symmetry 
breaking j2^,Q. In both the overlap and DWF 
case one can project out some number of low lying 
eigenmodes and correct the e{H{—M)) to ±1 |Q. 

The issue of whether there is a large number or 
even a finite density of small (zero) eigenvalues of 
H{—M) is a technical problem related to Wilson- 
like gauge actions since it can be shown that given 
bounds on a plaquette from unity there is also 
a corresponding bound on Xmin{H{—M)) |4^ . 
There are classes of configurations not satisfying 
these bounds that induce exponentially localized 
zero-modes of H{—M) which occur with non-zero 
density in a Wilson-like gauge action, and that 
have been observed to rapidly decrease with cou- 
pling The particularly troublesome modes 
are from dislocations of the gauge fields. The den- 
sity of these modes is greatly affected by changing 
the gauge action and lowering plaquette fluctua- 
tions [Q. In practice, the effects on e{H{—M)) 
can be fixed by projection or using a weaker cou- 
pling or both making the question of a finite den- 
sity somewhat moot. However, for a dynamical 
calculation in the chiral limit the zero level cross- 
ings of H{—M) should be repelled from crossing 
at M of the simulation, hence no ambiguity in the 
construction of a chiral fermion. 

At very strong coupling a new phenomena 
emerges where the density of the dislocations be- 
comes so large that the single flavor region in M 
merges with the doubler region (M > 2). Strong 
coupling and mean-field calculations differ over 
this mixing |4^ , and the approach to the contin- 
uum limit. 
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5. Conclusions 

Low lying fermionic modes provide a powerful 
probe of the vacuum. There are many studies 
using fermionic modes in quenched theories. Ob- 
viously more studies are needed with dynamical 
fermions. 

This work was supported by DOE contract DE- 
AC05-84ER40150 under which the Southeastern 
Universities Research Association (SURA) oper- 
ates the Thomas Jefferson National Accelerator 
Facility (TJNAF). 
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